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VO ! We review the main results of next-to-leading order QCD analyses of polarized deep- 

inelastic scattering data, with special attention to the assessment of theoretical uncer- 
tainties. 
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1. 

Experimental information on deeply inelastic scattering (DIS) of polarized leptons off 

polarized nucleon targets has greatly improved since the first E80/E130 experiments at 

SLAG (see ref. [Q for a complete bibliography.) While at an early stage the attention was 

Q^l mainly focused on the test of the Ellis- Jaffe sum rule, more recently the interest is con- 

^ ■ centrated on the study of the general features of polarized nucleons in the deep-inelastic 

region in the context of the QGD-improved parton model. The analysis of polarized DIS 

KA \ data can at present be performed using perturbative QGD at next-to-leading order accu- 

H \ racy, thanks to the recent computation [^ of order a| Altarelli-Parisi splitting functions in 

the polarized case. This analysis has been performed by many different authors |I[],0-|T^ 



with consistent results. I will present here the results obtained in ref. 0]. The strategy 
is the same as the one adopted in the case of unpolarized DIS: the polarized parton 
distributions Ag(a;, Q^), Ag{x, Q"^) at an initial scale Qo are assumed to have an arbitrar- 
ily chosen x dependence, specified by a set of unknown parameters; with the help of the 
QGD-improved parton model formulas and of Altarelli-Parisi evolution, one computes the 
structure function gi{x, Q"^) at each data point in the (x, Q^) plane, and fits the unknown 
parameters, which are in turn related to interesting physical quantities. 



A first important point is the test of the Bjorken sum rule |TT[]. The combination 
"1 



^^- /o 



dx 



g!{x,Q')-g-{x,Q') (1) 



can be shown to be proportional to the axial charge 

gA= [ dx [Am(x, Q^) - Ad{x, Q^)] , (2) 

Jo L J 

which is Q^-independent because of current conservation, times a Wilson coefficient 
Cns{c(s), which is perturbatively computable, and known to order a|. This is an im- 
portant and accurate theoretical prediction, since corrections to it may only come from 



isospin violation (~ 1%), from terms of order a| or higher in the Wilson coefficient, or 
from non-perturbative contributions, suppressed by powers of Ag^^/Q^. A direct test 
of the Bjorken sum rule has become possible, since gi data with deuteron and neutron 
targets are available. This is done by using the non-singlet axial charge qa as one of the 
parameters of the fitting procedure. The result |jl| is 

QA = 1.18 ± 0.05(exp) ± 0.07(th) = 1.18 ± 0.09, (3) 

to be compared with the value qa = 1.257 ± 0.003 measured in /5 decay (we will discuss 
in the next section the unceretainties in eq. (^. It can be concluded that polarized DIS 
data are consistent with the Bjorken sum rule at the level of one standard deviation, with 
an accuracy of less than 10 %. 

A second interesting question is the singlet contribution to the first moment of gi: 

1 
dxgi{x, Q^ 
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= d^\as)aom, (4) 

singlet 



where Cg (as) is the first moment of the singlet coefficient function, and ao(Q^) the 
singlet axial charge; oq is not scale independent because of the axial current anomaly. In 
the QCD-improved parton model, one can choose the factorization scheme 0] so that 
the first moment of the singlet combination of polarized quark densities, AE(1), is scale 
independent, and can therefore be interpreted as the total helicity carried by quarks. In 
this class of schemes one has 

aoiQ') = AS(1) - Uf^^Agil, Q% (5) 

where A(7(1,Q^) is the first moment of the polarized gluon density. The values of Oq, 
AS(1) and AglljQ"^) can be extracted from the fitting procedure outlined above. We 
find 

AS(1) = 0.46 ± 0.04 (exp) ± 0.08 (th) = 0.46 ± 0.09, 
A^(l,lGeV2) = 1.6 ±0.4 (exp) ±0.8 (th) = 1.6 ±0.9, (6) 

ao(oo) = 0.10 ±0.05 (exp) +0:1^ (th) = 0.10l°:i{. 

The results in eqs. (P) show that large values of AS(1) are compatible with small values 
of ao, provided Ag{l, Q^) is positive and large enough, as first suggested in refs. |jl2 . 



Finally, one can attempt using the value of as as one of the parameters in the fit |0] , as 
customary in unpolarized data analyses. It is interesting to note that the value obtained 
with polarized DIS data, namely 

asimz) = 0.120i°:°°t(exp)+0;00^(th) = 0.1201°:°^°, (7) 

is very close to other determinations, and that the uncertainty is reasonably small. 



We come now to a discussion of the theoretical uncertainties attached to the observables 
mentioned above, and summarized in table ^ This analysis has been presented in ref. [0], 
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AS 


A^ 


Oo 


Ots 


experimental 


±0.05 


±0.04 


±0.4 


±0.05 


+0.004 
-0.005 


fitting 

as & as 
thresholds 
higher orders 
higher twists 


±0.05 
±0.03 
±0.02 
±0.03 
±0.03 


±0.05 
±0.01 
±0.05 
±0.04 


±0.5 
±0.2 
±0.1 
±0.6 


±0.07 
±0.02 
±0.01 

+0.15 
-0.07 


±0.001 
±0.000 
±0.003 

+0.007 
-0.004 

±0.004 


theoretical 


±0.07 


±0.08 


±0.8 


+0.17 
-0.010 


+0.009 
-0.006 



Table 1 

Contributions to the errors in the determination of the quantities qa-, AS(1) 

A5f(l, IGeV^), ao(oo) and as{mz) from the fits described in the text. 
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o 
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where the interested reader can find more details. The experimental error is taken into ac- 
count by the minimum-square fitting procedure. We have added in quadrature systematic 
and statistic errors on each data point; this procedure results probably in an overestimate 
of the effective uncertainty on the fit parameters, since it does not account for correlations 
among systematics. 

A source of theoretical uncertainty which 
is often neglected is the arbitrariness in 
the choice of the functional form in x for 
the parton densities at the initial scale. 
We have considered a wide range of func- 
tional forms (see ref. [0] for details), and 
we have found the choice of the initial-scale 
parametrization affects consideralby the fi- 
nal results. In fact, different parametriza- 
tions lead to different estimates of contri- 
bution to the first moment of gi from the 
small-x region, where the experimental in- 
formation is very poor. The correspond- 
ing spread in the determination of physi- 
cal quantities must be included in the to- 
tal uncertainty. This is illustrated in fig. |l|, 
where the different curves refer to different 
parametrizations of the initial scale par- 
ton densities. The curves are quite close 
to each other in the measured region, as 

expected, since they all correspond to fits of comparable quality, while they differ consid- 
erably below X ~ 0.01. Perhaps, part of this uncertainty could be reduced using positivity 
constraints [T^ . 

Our analysis includes data points at Q^ down to 1 GeV^, in order to have a reasonable 
information at small values of x. At such low scales, one should worry about the uncer- 
tainty originated by higher orders in the perturbative expansion. These can be estimated 
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Figure 1: Plot of c/f(x,Q2) at Q^ = 10 GeV^ 
for different parton density parametrizations 
(labelled A-D). 



by varying the values of renormalization and factorization scales fiR, fip independently 
around Q^. Not surprisingly, this turns out to be the most important origin of theoretical 
uncertainty. It should be poited out that this "theoretical" uncertainty could eventually 
be removed, if more data at small x and higher Q^ were available; in fact, one could then 
exclude from the analysis data points below, say, Q^ ~ 4 — 5 GeV^, as in most unpolar- 
ized DIS analyses, thus avoiding the region where the QCD perturbative expansion is less 
reliable. 

Non-perturbative contributions are also potentially large at this low values of Q^, since 
they have the form of powers of Aq^j^/Q'^. Unfortunately, they are very difficult to 
estimate. One possible strategy is that of comparing results obtained fitting all data 
above Q"^ = 1 GeV^ with those obtained by excluding data points below 2 GeV^. This 
procedure indicates that the contribution of power-suppressed terms is not very large, 
compared to other sources of uncertainty. A similar conclusion is obtained by fitting the 
Bjorken sum to its perturbative expression, supplemented with a twist-4 term a/Q"^, with 
the parameter a taken from renormalon and sum rule estimates. 

Other minor sources uncertainties, such as violations of the SU{3) fiavour symmetry or 
the position of heavy quark thresholds in Q^ evolution, are also listed in table |I]. 

I wish to thank G. Altarelli, R. Ball and S. Forte for the fruitful collaboration on the 
subject of this talk. 
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